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Abstract 

C^*) ■ In this article, we assume that the two quarks have unequal masses and calculate the next- 

to-leading order contributions to the spectral densities of the two-point correlation functions of 
the vector, axial-vector, scalar and pseudoscalar currents. We take dimensional regularization 
to regularize both the ultraviolet and infrared divergences, and use optical theorem to obtain 
the spectral densities directly, furthermore, we present some necessary technical details for 
readers convenience. The analytical expressions are applicable in many phenomenological 
analysis besides the QCD sum rules. 
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Q-e 1 Introduction 

p ,1 

The QCD sum rules is a powerful theoretical tool in studying the ground state mesons [TJH], and 
has given many successful descriptions of the meson properties. In the case of the interpolating 
currents consist of quarks with equal masses or zero masses, the leading order and next-to-leading 
' order perturbative contributions to the spectral densities of the correlation functions are known 

^ , [2J [3] . There are few works on the next-to- leading order perturbative contributions to the spec- 

' tral densities of the interpolating currents for unequal quark masses, as the calculations are very 

difficult. 

In Refs.[H[5], Reinders, Rubinstein and Yazaki calculate the next-to- leading order contributions 
to the spectral densities for the pseudoscalar and vector currents in the unequal mass case. Some 
technical details are presented in Ref.[3]. In the famous review (see Ref.[3]), Reinders, Rubinstein 
and Yazaki illustrate how to use optical theorem to calculate the next-to-leading order contributions 
to spectral densities of the correlation functions (for the scalar, pseudoscalar, vector and axial- 
vector currents) directly, and regularize the divergences in the massive-gluon scheme. In Ref.[6], 
Schilcher, Tran and Nasrallah also resort to the massive-gluon scheme to regularize the infrared 
divergence, and use optical theorem to calculate the next-to-leading order contributions to the 
longitudinal and transverse spectral densities directly for arbitrary conserved or non-conserved 
vector or axial- vector currents of the massive quarks, and obtain analytical expressions. Although 
those works are all performed in the massive-gluon scheme, the resulting (infrared) divergent terms 
are different from each other, the expressions also different from each other. 

In Ref.[7], Generalis calculates the two-loop Feynman diagrams directly using the dimensional 
regularization, and obtain the next-to-leading order contributions to the correlation functions for 
the vector and axial-vector currents for unequal quark masses. We have to obtain the spectral 
densities through dispersion relation from the cumbersome expressions, it is a difficult work. 

In this article, we resort to optical theorem (or Cutkosky's rule) to calculate the next-to-leading 
order contributions to the spectral densities of the correlation functions for the scalar, pseudoscalar, 
vector and axial- vector currents in the case of unequal quark masses. In calculations, we regular- 
ize both the ultraviolet and infrared divergences in the scheme of dimensional regularization, and 
present the necessary technical details for the readers convenience, i.e. one can check the calcula- 
tions. There are two routines in application of optical theorem, we resort to the routine used in 
Refs.[U|B], not the one used in Rcf.|8j. After the present work is finished and submitted to the net 
http: / / arxiv.org/[ some friends kindly draw my attention to Refs.[51 ITU]. In Ref.[5], Djouadi and 



Gambino study the QCD corrections to the electro-weak gauge boson self-energies and the Higgs 
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self-energies for arbitrary momentum transfer and for different internal quark masses, and obtain 
explicit expressions for both the real and imaginary parts of the self-energies. They calculate the 
two-loop Feynman diagrams directly using the dimensional regularization, the tedious work is dif- 
ficult to follow. On the other hand, if we calculate the imaginary parts of the correlation functions 
(or the self-energies) directly via optical theorem, the calculations are greatly simplified and more 
easy to follow, furthermore, we can separate the contributions come from the real and virtual gluon 
emissions explicitly, which have direct applications in studying the decays of a polarized W boson 
into massive quark-antiquark pairs, such as t — > bW + — > bqiq 2 , H — > W~W* + ~ > \V~q1q2, the H 
denotes the Higgs boson [ID] . 

The article is arranged as follows: we calculate the next-to-leading order contributions to 
spectral densities in Sect. 2; in Sect. 3, we present the final analytical expressions; and Sect. 4 is 
reserved for our conclusions. 



2 Explicit calculations of the spectral densities at the next- 
to-leading order 

In the following, we write down the two-point correlation functions in the QCD sum rules, 

n^(p) = i /dW^<o|T{j M (x)4(o)}|o), 



P tl p V \ TT \P^P V 



n<j/p(p) = — — 2 n s/P (p) , 

Il s/P (p) = * /Vxe*-*(0|T{J(x)jt(o)}|0), (1) 

where J^(x) — Jy(x), J^(x) and J{x) — Js(x), Jp{x). The lower indexes denote the scalar (5), 
pseudoscalar (P), vector (V) and axial- vector (A) currents respectively. The correlation functions 
can be expressed in the following form through dispersion relation, 

TT Ja 2 S- P 

where i — S,P, V, A, the A 2 = (mi + m,2) 2 is the threshold parameter, and 

= Pi (s) + P t (s) + pi (s) + • • • , (3) 

the Pi{s), pj(s), Pi(s), ■ ■ ■ are the spectral densities of the leading order, next-to- leading order, 
and next-to-next-to-leading order, • • • . At the leading order, the perturbative spectral densities 
are rather simple, 



3 VA(s,m 2 ,m 2 ) J 2 A(s,m 2 ,m^) 
Pv ' a{s) = 8^ a r (WlTm2) 3s 



where 



n / x (wi =F mo) 2 3 v A(s, m?, mf) r . , 9 , ... 

P°S/p(*) = , ^ V / 2; { S -(mi±m 2 ) 2 } , (4) 



W 2 2 \ 2|4|4r>2r>2o^2 / r \ 

\{s,m 11 m 2 ) = s + m 1 + m 2 — Zsm 1 — Zsm 2 — Zm 1 m 2 . (5) 



At the next-to-leading order, there are three Feynman diagrams contribute to the correla- 
tion functions, see Fig.l. We calculate the imaginary parts (or the spectral densities) using the 
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Figure 1: The next-to- leading order contributions to the correlation functions. 




Figure 2: Six possible cuts correspond to virtual gluon emissions. 



Cutkosky's rule or optical theorem, the two approaches lead to the same results, then use the dis- 
persion relation to obtain the correlation functions. There are ten possible cuts, six cuts attribute 
to virtual gluon emissions, see Fig. 2, and four cuts attribute to real gluon emissions, see Fig. 3. 

2.1 Contributions from the virtual gluon emissions 

The six cuts shown in Fig. 2 attribute to virtual gluon emissions and correspond to the self-energy 
corrections and vertex corrections respectively. The first four cuts corresponding to one-loop 
fermion's self-energy corrections, we calculate the Feynman diagrams directly using the dimensional 
regularization and choose the on-shell renormalization scheme to subtract the divergences so as to 
implement the wave-function renormalization and mass renormalization. Then we can take into 
account all contributions come from the six cuts shown in Fig. 2 by the following simple replacement 




Figure 3: Four possible cuts correspond to real gluon emissions. 
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for each vertex in the interpolating current^, 

u(pi)-l^u(p 2 ) -s- u(p 1 )j fl u(p 2 ) + u(p 1 )T fl u(p 2 } 



Z 1 y / Z 2 u(pi)^ pi u(p2) +u(pi)T fi u(p 2 ) 

u(pih^u(p 2 ) (l + ^5Zi + ^SZ 2 j + G(pi)r ^u{p 2 ) , (6) 



where 



% = ^ + ^ = l + i^(--^ + -^ + ho g ^t + 3 - 1 -l) , (7) 
3 7T V 4e uv 2e IR 4 & Anfj, 2 4 / w 

are the wave-function renormalization constants come from the self-energy corrections, see Fig. 4, 
and 

r 4 2(' /• 1 - x , fd D k f 

T M = 7mo3s / dx I d v 



r s J J »J (27T) 

r ( 3 ) J 7,2 / , \ , ,.„ 1 1 -1 , 



[k% + {xpi + yp 2 ) 2 ^ \ a \ - 



*b 1 - n £ uv + 2(1 - z)(l - t/)(s -mi- m 2 ) 



4 /•! 



+2{x + y)mxm 2 + 2x(l - x)m\ + 2y(l-y)ml} + -g 2 s i dx I dy 

3 JO JO J (^TTJ 

4r(3) { [x 2 mi - 2/(1 - x)m 2 ] pi M + [y 2 m 2 - z(l - y)m 1 ] p 2f _ l } 
[k% + {xpi+yp 2 ) 2 } 3 



(8) 



comes from the vertex corrections, see Fig. 5. Here 7 is the Euler constant, /x 2 is the energy 
scale, and the Euclidean momentum fog = (fci, fc 2 , £3, /C4). In this article, we take the dimension 
D = 4 — 2euv = 4 + 2eir to regularize the ultraviolet and infrared divergences respectively, and 
add the energy scale factors /i 2euv or /i _2eiR when necessary. In the limit mi = m 2 = m, the T M 
can be reduced as 

4 2 r 1 z -1 -* 

F M = 7mo#s / da; / d v 



r s J Q J "J (2n) D 

m ' ' fe| fl - ^ uv ) + 2(1 - (x + y) + xy)(s - 2m 2 ) 



[k% + (xpi + yp 2 ) 



4 rl rl-x 



+4(x + y)m 2 - 2(x + y) 2 m 2 + Axym 2 } + -g 2 / / ,. 

2r(3)m [{x + y) 2 - (x + y)] (p lM +p 2A1 ) 
[&I + (xpi + yp 2 ) 2 ] 3 



(9) 



then the resulting expression is greatly simplified. 

We carry out the integral over the variables x, y and fee, and observe that the ultraviolet 
divergent terms in the T^, 8Z\ and <5Z 2 are canceled out with each other, the net contributions 
have no ultraviolet divergence, 



Fp = \-l»f{ S ) + \-h{s)pi» + \-f2(s)p 2 », (10) 



"Here we use the vector current to illustrate the procedure. Other currents can be treated analogously. 
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Figure 4: The quark self-energy correction. 



where 

„, , . 2 , 47ru 2 , mim 2 2(s — m? — m%\ , 

f(s) = f(s + — +log— ^ + 2 7 -3 + 31og- 1 -^- 1 1 2 / log 
£ir s 4tt/^ ^A(s,m 2 ,m|) 

1 

7±(«) = F(s)+2( S -m 2 -m 2 ) [Foo(s)-^io(s)-^,i(s) + yii(s)] ±2 mi m 2 

[V 10 («) + V 01 (a)] + 2m 2 [V 10 (s) - F 20 (s)] + 2m 2 [V i(s) - V 02 (s)] , 
/i±(s) = 4miV2o(s) T 4m 2 Vbi(s) ± 4m 2 Vn(s) , 
/ 2 ±(s) = ±4m 2 Fo2(s) - 4mxVxo(s) +4mxFxx(s) , 

ys— (mi — m 2 j^ 

and s = p 2 , the V(s), Voo(s) and Vy(s) with i,j = 0, 1, 2 are given explicitly in the appendix. Here 
we introduce the lower-indexes ± to denote the corresponding expressions of the vector (+) and 
axial- vector (— ) currents respectively. The infrared divergence ^ comes from the wave-function 

renormalization constants Zi in Eq.(7), while the infrared divergence log ^yz^j) ^ comes from 
the term 

1 rl-x 



4 2 f , [ . /" d C *JS 2r(3)(.s -m 2 -m 2 ) 



7m - 5 2 dx dy 7^—^ ^— § (12) 

in Eq.(8). The infrared divergent terms obtained in the present work are consistent with that 
obtained by Groote, Korner and Tuvike in Ref.[10]. In the massive-gluon scheme, we expect that 

( \ m 2 m 2 

the infrared divergences — and log ( ] — appear in the forms log -, r 2 — yy (or log — — 1 and 

° £ir ° V 1— w j £ir L L ° (mi+m2j v ° raim2 f 

log log ( m "+n l2 )i ( or lo § (t=£j) 1o S wfe)> respectively, where the m ff is the supposed gluon 

mass. In Ref.[2], there only exists the infrared divergence log while in Ref.[6], there exist 

(\ 2 2 

) log mg and log — log ™ a . In this article, we calculate the 
l—UJ J & mim2 & m2 & mima 1 

vertex corrections using the conventional Feynman parameters, while in Ref.jB], Schilcher, Tran 
and Nasrallah calculate the vertex corrections using optical theorem and dispersion relation. 

The total contributions come from the virtual gluon emissions (see Fig. 2) to imaginary parts 
of the correlation functions can be expressed in the following form, 

Imn^( S ) 4a s 3 [ dP^px d D ^p 2 D D 

- ' - (27r) 5 (p-pi ~P2) 



TV 



3 7T ir J {2tt) d - 1 2E Pi {2tt) d - 1 2E, 



l>2 



s — (m 1 — m 2 ) 2 
P\nP2 V + VipViv 2 3liu 



fi(s)(mipifj,P2u - rn2PifiPiv) 



-f2(s)(m 2 p 2 f_ l pi )} ■ (13) 
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Figure 5: The vertex correction. 



In this article, the up-indexes V and R denote the contributions come from the virtual and real 
gluon emissions respectively. We carry out the integrals in Eq.(13) directly in D = 4 + 2eiR 
dimension as there is no ultraviolet divergence, and obtain the following results, 



ImIl£ M ( S ) 



4a, 

3 7T 



Pv/a(s) I ^- - 2 log47r + 2 7 - - + - lo ; : 



7 1 \ 2 (s,m\,m\)m\ml 



/J, 8 S 3 



s — m? 



n 2 ) 



log 



1 -OJ 



+ 4a s y/\{s,m\,ml) 



1 



3 7T s 2 1 12tt 2 

(mi ± m 2 ) [/i±(s) + /2±(g)] 
32tt 2 



1 

£lR 

1 - 



2 log 4?r + 27 - 2 + log 



A(s, m 2 , to 2 . 



2\1 



s — m? — m. 



loe 



1 + w 

1 -UJ 



(14) 



Imn^ /p (s) 



4 a 
3^ 

s 



S P°s/p( s ) 



f 1 7 

J 2 log 4^ + 2 7 - - 

Uir 2 



1 i„„ A 2 (5,m 2 ,TO 2 )mfm3 1 

2 l0S ^ + 2 /±(5) 



dog 



1 



1 -U> 



1 „ A(s,m?,m 2 ) 

2 log47r + 2 7 - 2 + log v 21 

£ir fJr 



^(s,™ 2 ,™ 2 ) 

4 a s 3 y A(s, to 2 , m 2 ) r 2 2w 2 2\ r/ i \ \ c 1 \ 1 

+ q — 2^ U s + m i ~ m 2)( s + m 2 - TO i) [/i±(s)mi ± f 2 ±(s)m 2 ] 

Tfi±{s)m 2 {s + ml-ml) 2 - f 2 ±{s)m 1 (s + ml-m 2 1 ) 2 } , (15) 

where we have used the 75 symmetry of the interpolating currents to obtain the spectral densities 
^n v A (s) and imn£( s ) _ 



2.2 Contributions from the real gluon emissions 

The four cuts shown in Fig. 3 correspond to real gluon emissions. The scattering amplitudes for the 
real gluon emissions are shown explicitly in Fig. 6. From Fig. 6, we can write down the scattering 
amplitude T° Q (p), 



2 ft — mi 



7/j > I'p. 



~i>i- fi-m 2 



Y la j v 



(P2) , (16) 
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where A a is the Gell-Mann matrix. Then we obtain the corresponding contributions Imlly(s) and 
Imllg (s) to the imaginary parts of the correlation function with optical theorem, 



Imll£(s) 



1 r d D ~ 1 k 

D- 1) J (2n) D - 1 2i 



jD-1 



Pi 



d 



D-l 



P2 



2tt(D -I) J {2n) D - 1 2E k (2w) D - 1 2E pl (2tt) d - 1 2E P2 



(2w) D S D {p-k-p 1 -p 2 ) 



tt(D-I) 



jD-1i 



P 

d D ~\ 



lD-l 



P2 



(2ir) D -^2E k {2n) D - 1 2E pl (2ir) D ^2E P2 

s-K 2 1 



{2Tt) u 5 u {p-k- Pl - P2 ) 



m\ 



s — m\ — m\ 



{k-pi) 2 {k-p 2 ) 2 
X(s,m 2 ,m 2 ) 



+ 



k-p\k-p 2 k-p 1 k-p 2 



3s 

{s-K 2 ) 2 
k -pik ■ p 2 



2 + 



+ 4eiR [s - (mi - m 2 ) 2 ] 
(mi - m 2 ) 2 ' 



mr 



[s - (mi 

™2 



m 2 ) 



{k-pi) 2 {k-p 2 ) 2 k-p\k-p 2 



+ 16 



(17) 



Imnf( S ) 



1 

2^ 



J2_gi_ f d"- 1 
it J {2w) D - 1 

{ 2[ 



d D - x k 


d°- 


¥i d D ~ 


l P2 


(2n)D-i2E k 


(2ir) D - 


12E P1 (2ir) D - 


-12E P2 


(pKl( P )}g 


a ^P V 
p 2 






• d D ^k 


d D 


~ l Pi d D 


-'P2 



(2ir) n 5 n (p-k- Pl -p 2 ) 



(mi - m 2 f 
s 

s-K 2 
k-pik-p 2 



2E k {2tt) d - 1 2E p1 (2k) d ^2E P2 

2 

(mi+m 2 ) 2 ] - " n - i 



{2n) D 5 D {p- 

™2 



k-pi -p 2 ) 



m\ 



m% 



(k-pi) 2 (k-p 2 ) 



k-p 1 k-p 2 



(s - K 2 ) 2 
k-p x k-p 2 



(18) 



where we have used the identities ^2 u{p\)u(p\) =fe\ + mi and ^v{p 2 )v{p 2 ) =t>2 — m 2 for the 
particle and antiparticle respectively, and take the notation K 2 = (pi +p 2 ) 2 . We carry out the 
integrals in Eqs. (17-18) in D = 4 + 2sm dimension as there is only infrared divergence, and obtain 
the contributions come from the real gluon emissions (see Fig. 3) to spectral densities, 



i m n« M ( s ) 



4 a 
3^- 



S Pv/a(s) \ -— + 21og47r - 2 7 + r _ bg 



A 3 (s, m\,m 2 ) 



m\m\s 2 pfi 



+ (s - m\ - m\)R 12 (s) 



-Ru(s) - R 22 (s) - 



2(s-m\-m 2 ) fl + uj 



3\/A(s, m 2 , m\) 



log 



R{ 2 (s) + 



log 



1 



21og47r + 27-2 + log 

£ir 



1-w 

A 3 (s,m 2 ,m|) 



s — mr — mS 



m\m\s 2 [i^ 



4a s 

3 7T 



,s - (mi =F m 2 ) 2 



4S7T 2 



log ( Y~ ) ( s ~ m i ~ m 2) _ \/A(s,m 2 ,m 2 ) 



+ 



16s7T 2 



i? 2 2 ( S ) 



2 + 



(mi =F w 2 ) 2 



(19) 
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Imnf /p (s) 4a 



Figure 6: The amplitudes for the real gluon emissions. 

\ 3 (s,ml,m 2 ) 



i— Ps/p( s ) ^ - — + 21og47T - 2 7 + 2 - log " ^ff¥- +(s-mi- m 2 2 )R 12 (s) 
3 7r ' I £ir m^m^s /z 



-Ru(s) - R 22 (s) - R{ 2 (s) + 



R 



12 



+ 



mt — m. 



log 



1 + U) 
1-uj 



2 s-(mi±m 2 ) 2 V /A(s,mf,m^) 
A 3 (s,mf,m|)" 1 



21og4^ + 2 7 -2 + log ., ., ., 

£ir m^m^s^/U 4 



(20) 



the expressions of the i?n(s), i?22(s), -Ri2(s), -Ri2(s) and R\ 2 (s) are given explicitly in the ap- 
pendix. We have used the 75 symmetry of the interpolating currents to obtain the spectral densities 
imn%(s) and Imng( S ) ^ 



3 Analytical expressions of the total contributions 

We add the contributions come from the virtual and real gluon emissions together, and obtain the 
total perturbative spectral densities at the next-to-leading order, 



4 a 



Pv/a( s ) = ^—Pv/A( s )\^f±( s )- R ii( s )- R 22(s) + {s-m 1 - m2 )R 12 (s)- 



+2 log ■ 









\f± 




m 2 s 


2 2\ 


,mi,mi) 



2(s — m\ — to 2 ) 



log 



2(s — mj — m 2 ) 

l+LO 



log 



LO \ , A(s, TO?, TO9) 

log 



TO1TO2S 



+ 



+ 



3- v /A(s,m?,m|) V 1 
4a s f s — (mi =F m 2 ) 2 



^i 2 («) 



3 7T 



4S7T 2 



^/A^TO 2 ,??! 2 ) 
16S7T 2 

1 /_ s- 



R 2 12 (s) 



log ( s - m i - m l) - \Jx{s,mj,ml) 

(mi Tm 2 ) 2 



2 + 



1 



log 



12^ 2 I ^(s,™ 2 ,™ 2 ) VI 



s 

1 + u> 



(mi ± TO 2 ) (/i±(s) + f2±(s)) 
32tt 2 



(21) 
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p\ /P {a) = ~Ps /P (s)|^/±(s)- Rn(s)- R 2 2(s) + (s - ml - m 2 2 )R 12 (s) - ^ 

\J ml mis 2(s — m? — to 2 .) , + A(s, m?,mo) 

+2 log ~ /_ - ^ 2 2 ^ + -^== = o 2 / log ( -)log v 



\(s,m(,m%) V /A(s,m 2 ,m|) V 1 -^/ ™i m 2S 

pi r ^ i ^(*) 1 1 

-«i 2 W + 2 S -(m 1 ±m 2 ) 2 / 

+ o — oo 2 2 {(s + m 2 - to 2 )(s + to 2 - m 2 ) [/i±(s)mi ± / 2 ±(«)m 2 

T/i±(s)m 2 (s + m 2 -m^) 2 -/ 2± (s)mi(s + m^-m 2 ) 2 } . (22) 



The infrared divergences log ^y^j J ^ and constant terms log47r, 7, log^i 2 come from the 

virtual and real gluon emissions are canceled out with each other. The spectral densities at the 
next-to-leading order do not have explicit energy scale dependence, although the strong coupling 
constant a s (/j.) is scale dependent quantity. In the appendix, we present some necessary technical 
details for readers convenience, one can follow the necessary steps to check the calculations if one 
does not feel confident with the final expressions. 



4 Conclusion 

In this article, we calculate the next-to-leading order perturbative contributions to the spectral 
densities of the two-point correlation functions of the vector, axial-vector, scalar and pseudoscalar 
currents. In calculations, we assume that the quarks have unequal masses and use optical theorem 
to obtain the spectral densities directly, furthermore, we take the scheme of the dimensional regu- 
larization to regularize both the ultraviolet and infrared divergences. The ultraviolet and infrared 
divergent terms are canceled out with each other separately, the net spectral densities are free of 
divergences. The analytical expressions are applicable in many phenomenological analysis besides 
the QCD sum rules. 
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Appendix 

Firstly, we write down the fundamental integrals in calculating the vertex corrections, 



v ab{s) = iw^rv^fe * vr(3) 3 , 

{) Jo Jo " J (2n)V [k 2 E + {xpi+yp2) 2f> 



E 

2 



fc|r(3) 
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and carry out the integrals to obtain the following analytical expressions, 

^oo(a) = ,„ \ I- ^ l°g T^+lV 10 ^"^ - log 2 d + <*) 



y/\(s,ml,ml) I \l-wy\£i R 4tt/x 2 
+ l0g2(1 r" 22) - log 2 (l + c 2 ) + 2 log( Wl + W2 ) log f l±HL) toga* log f i±^) 

- lo ^ 2 10 § (l^) - Li2 (iT^) Li2 ( TT^) + "'I ' 



= Voo(s) - t= L= log ( ) ( — + log -i-j + 7 



V W {s) = -{-log 1 -—log +log^ 

s ^ 2 \ 1 — J u 2 \ 1 — to J cji 

Vbi(s) = ^io(s)Liow 2 . 

V20O) = tt 1 ; — lo s 1 ? — ; — v lo s 1 H ; — lo s 1 2 

+ — r~ log — + 1 \ , 

Vo2(s) = V2o(s)| Wl+ + W2 i 

^u( s ) = 7T1 — : — lo s 1 — : — x i°s n 2 ; — lo s h 

2s [LJ1+LJ2 \1 - w / 2(cji+u;2) \l _w 2/ W1+W2 
wi wi w 2 w 2 1 

h ■ — log — 1 ■ — log n , 

UJ\ + u> 2 W 2 OJl +0J2 0J\ J 

V s = +log— ^-7+1 log logl-w 2 

euv s LO1+LO2 \l -uj J wi +lu 2 

w l , 2n Wll0gWi + W 2 logCJ 2 . 01 / . s 

log(l - LO2) - 2 ■ h 2 log(wi + w 2 ) , 



LUl +UJ 2 U)\ + UJ2 

= V(s) + — +log^- 7 + l, (24) 



where 



Wi = 

w 2 = 
M = 



s + m 2 — m\ 

s + m 2 — to 2 
mi + TO2 



mi — m2 



Jo 



U,{.r) = - / dt^lzl). (25) 



Secondly, we take the notation 

iD-lZ jD-1~ jD-l; 



(a f d^kd^pid^p2 D 

J dps = 7 2^ 2 ^ 2 ^ 2 J (p-*-*-*). 
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for simplicity, and write down the analytical expressions of the three-body phase-space integrals, 



R (s) 



dps 



TO? - TOn , / M + LJ 

— —. — - lo™ ' 



\M -LJ 



sm\ + srrvk — 2m\mk ( 1 + lj 
log 



4s 



+ 



yj\{s,m\,m\)(s + m\ + m|) 



8s 



1-LJ 
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where s = s — (mi — m 2 ) 2 . 

In the following, we present some necessary technical details in calculating the integrals for 
both the virtual and real gluon emissions, one can check the calculations; here we smear the lower 
index IR in ein for simplicity. 
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where 



ui = 



u 2 
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K 2 + m,2 — m\ ' 
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and we have neglected the imaginary terms. The infrared divergent terms of the form log ^j^j^ \ 

come from the diagrams of the virtual gluon emissions (vertex corrections) and real gluon emissions 
are canceled out with each other. 

In calculating the divergent integrals, we have used the following trick, 
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the functions f(x,e) and g(;r) have no poles in the range x = — lj and e — > + . 



(w - a;) 1 - 5 

(30) 
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